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Abstract 

In the absence of coalescence, coarsening of emulsions (and foams) 
is controlled by molecular diffusion of the dispersed phase species from 
one emulsion droplet (or foam bubble) to another. Previous studies 
of dilute emulsions have shown how the osmotic pressure of a trapped 
species within droplets can overcome the Laplace pressure differences 
that drive coarsening, and "osmotically stabilise" the emulsion. Web- 
ster and Cates (Langmuir, 1998, 14, 2068 — 2079) gave rigorous crite- 
ria for osmotic stabilisation of mono- and polydisperse emulsions, in 
the dilute regime. We consider here whether analogous criteria exist 
for the osmotic stabilisation of mono- and polydisperse concentrated 
emulsions and foams. We argue that in such systems the pressure dif- 
ferences driving coarsening are small compared to the mean Laplace 
pressure. This is confirmed for a monodisperse 2D model, for which 
an exact calculation gives the pressure in bubble i as Pi = P + U + Pf, 
with P the atmospheric and IT the osmotic pressure, and Pf a 'geo- 
metric pressure' that reduces to the Laplace pressure only for a spher- 
ical bubble, and depends much less strongly on bubble deformation 
than the Laplace pressure itself. In fact, for Princen's 2D emulsion 
model, Pj is only 5% larger in the dry limit than the dilute limit. We 
conclude that osmotic stabilisation of dense systems typically requires 
a pressure of trapped molecules in each droplet that is comparable to 
the Laplace pressures the same droplets would have if they were spher- 
ical, as opposed to the much larger Laplace pressures actually present 
in the system. We study the coarsening of foams and dense emulsions 
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when there is insufficient of the trapped species present. Various rate- 
limiting mechanisms are considered, and their domain of applicability 
and associated droplet growth rates discussed. In a concentrated foam 
or emulsion, a finite yield threshold for droplet rearrangement among 
stable droplets may be enough to prevent coarsening of the remainder. 



1 Introduction 



Previous work on dilute emulsions has shown that, when coarsening is solely 
caused by a diffusive flux of dissolved dispersed-phase molecules between droplets 
(Ostwald ripening |K|), coarsening may be prevented by the addition of a suffi- 
cient number of molecules that are insoluble in the continuous phase and hence 
trapped within droplets. The trapped molecules provide an osmotic pressure 
which counteracts the Laplace pressure due to surface tension [[J (which drives 
coarsening); resulting in 'osmotic stabilisation'. A quantitative criterion for sta- 
bility, valid even for emulsions with polydispersity in both the droplets' sizes and 
number of trapped molecules they contain, was given in |9| . It was also found in [[J 
that for an 'insufficiently stabilised' emulsion (without enough trapped species to 
obey the required criterion), the subsequent coarsening was qualitatively unal- 
tered from that without any of the trapped species [fTl, |12]l . The only effect was is 



to reduce the effective volume fraction of the coarsening droplets by an amount 
corresponding to the final volume of a population of small droplets that attain co- 
existence with the coarsened bulk phase. The latter are prevented from entirely 
dissolving by the trapped molecules they contain. When the effective volume 
fraction is reduced to zero, full stability is achieved. 

In the present paper we discuss whether analogous conditions exist for the osmotic 
stabilisation of concentrated emulsions and foams containing trapped insoluble 
molecules in the dispersed phase; we again consider both the monodisperse and 
the polydisperse case. Previously the use of Nitrogen as a trapped species to 
stabilise foams has been investigated with theory |14|, [Tj| , experiment p4|HT? 



and computer simulations [0, but theoretical conditions for stability of foams 
with nonspherical bubbles were not found. The use of other trapped gases to 
inhibit dissolution (as opposed to coarsening) of spherical bubbles is addressed 
in |TJ| where the effect of condensation of the included gases is also addressed; 
we do not consider this here. In our work, we treat the idealised case of a fully 
insoluble trapped species (perfect trapping) whereas in many cases the results will 
be modified by residual solubility effects (considered for the dilute case in 0). 
Some gases are, however, practically insoluble in water, for example C2F 6 [pOfl ; and 



in the case of emulsions, the insoluble limit is easily achieved (by using oligomeric 
species) . For simplicity we do not consider the effects of residual solubility in the 
present work. 
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From now on we use mainly the language of foams, in which the trapped and 
soluble gas are treated as ideal, but the work is also applicable to dense emulsions 
whose droplets contain an ideal mixture of soluble and trapped molecules. The 
only distinction between these cases is that foams are compressible, but in fact 
much of the paper concerns effectively incompressible foams, in which pressure 
differences between bubbles are negligible compared with their mean pressures, 
and the total bubble volume is conserved. 

The osmotic stabilisation of foams is a far more complex proposition than for 
dilute emulsions. In a dilute emulsion, each droplet is spherical with an internal 
pressure that is increased by the droplets Laplace pressure, which is directly 
related to the bubble's volume through its spherical geometry, and equals 2a/ Ri 
for a droplet of radius Ri and surface tension a. Since nonspherical foam bubbles 
in contact have no simple relationship between their volume and surface area, 
there is no direct relationship, in a concentrated foam, between a bubble's pressure 
and its volume. 

A bubble i with a radius of curvature at its Plateau borders has a Laplace 
pressure of 2a/ri, an d an internal pressure of Pi = P + 2a /ri. Given a connected 
domain of the liquid phase, the atmospheric pressure P equals the pressure in 
the Plateau borders (we ignore the effects of gravity throughout). The drier 
the foam, the smaller the borders and the larger the Laplace pressure, so that 
if osmotic stabilisation of a foam or emulsion required the pressure of trapped 
molecules to balance the Laplace pressure (as it does in the dilute case) osmotic 
stabilisation would be hard to achieve. However, another relevant length scale 
is Ri = (3\^/47r) 1 / 3 , where Vi is the volume of bubble i. Since the curvature of 
adjacent bubble-bubble faces are of order 1/Ri, the pressure differences between 
adjacent bubbles (which are responsible for coarsening), are of order a/Ri <C cr/Vj. 
So if the osmotic stabilisation of a foam merely requires the partial pressure of 
trapped molecules to balance a/Ri, then osmotic stabilisation is a reasonable 
proposition. We shall confirm below that this is so. The reduction in the driving 
force for coarsening from the level suggested by the Laplace pressure has long 
been recognised (see e.g., [T3[]) but the implications of this for stabilisation by 



trapped species has not previously been addressed in detail. 

In what follows we first give in Section ^] some results for the dilute case (based 
on 0), and then investigate what factors determine the pressure within concen- 
trated foam bubbles. We follow the approach of Princen [pi|-p3|, and study the 
osmotic compression (at atmospheric pressure P) of previously spherical foam 
bubbles by an osmotic pressure II. Section |3] discusses the disjoining pressures 
between bubbles, the uniformity or otherwise of the osmotic pressure II, and the 
condition for mechanical equilibrium with an excess bulk gas phase created by 
coarsening. The increase in a bubble's pressure above that of a bulk gas is defined 
by Pj G = Pi — P — LT, and we argue that typically PP ~ a/Ri. This is explicitly 
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confirmed in Section |4] for Princen's monodisperse 2D model. We call Pf the 
geometric pressure and identify it (rather than Laplace pressure) as the driving 
force for coarsening, in general. A formal condition for osmotic stabilisation is 
derived in Section |5], and examined for various limiting geometries. We discuss 
Pf for polydisperse foams in Section |], and consider exceptions to our estimate 
that Pf ~ c/Ri which might arise under certain conditions (which we argue to 
be uncommon). We then use the geometric pressure to further investigate the 
stability requirements of polydisperse foams. 



The coarsening of insufficiently stabilised foams is studied in Sections 0-|TT|, start- 
ing with a simple mean-field model in Section [8[ By considering dissipation rates 
for the diffusive flux of gas between bubbles and for the viscous rearrangement 
of bubbles, we are able to predict (as a function of various parameters) the rate- 
limiting mechanism and the associated growth rate (Section ||). Section |10| shows 
that when bubble rearrangements are sufficiently rare, elastic stresses may arrest 



coarsening. Section |TT] extends this to the case of a finite yield strain beyond 
which bubble rearrangements will cause the foam to flow, finding that a foam's 
initial state then determines whether coarsening will still occur. We conclude in 
Section with a brief summary and discussion of our results. 



2 Dilute Foams (Spherical Bubbles) 

A "dilute foam" comprises of spherical gas bubbles floating freely in a solvent. 
Following 0, we treat the bubbles as macroscopic objects and neglect their en- 
tropy of translation. 

Firstly we consider the size (and hence composition) at which spherical bubbles, 
containing both soluble and trapped gas molecules, may coexist with a bulk gas 
phase (for example, formed by one bubble in the foam becoming macroscopic). 
The gas pressure within a spherical bubble labelled i is increased above the 
atmospheric pressure by its Laplace pressure 2a /Ri\ 

p. = p + 2 i a) 

rii 

We consider bubbles with Nf soluble gas molecules and Nf trapped gas molecules, 
and treat the gases as ideal. Then Pj = Pf + Pf, where Pf, Pf are the partial 
pressures of the soluble and trapped gas molecules respectively, and for the soluble 
gas molecules 

ti = ri + kTln(%) (2) 



P 

where \l\ is the chemical potential of a bulk gas phase of soluble molecules at 
atmospheric pressure P. Using Eqs.[l]j2| and Pf = p — P/\ we may write A/ij = 
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Hi — nl as 

A^ = fcTln(l + 2(T/ ^^ T ) (3) 

So that when Pj = 2a / R iy /^f = [i s b and bubbles may coexist with a bulk gas phase 
(at pressure P). Using the ideal gas law for the trapped species PjVi = NjkT, 
such coexistence requires 

2a NfkT 

Ri (4tt/3K 3 1 1 

where Vi = (47r/3)i?f. This expression determines a 'coexistence volume', , 
at which the Laplace pressure and the partial pressure of trapped gas balance. 
Solving for V { B we have 

This is identical to the coexistence volume for a dilute incompressible emulsion 
droplet, of interfacial tension a, containing Nf molecules of a trapped species ||. 

At coexistence with a bulk gas phase, A//j = and P? = P 6 S = P, so from the 
ideal gas law the number of soluble gas molecules in bubbles is 

kT kT\Aix 1 2a 



Similarly in two dimensions, dilute (circular) bubbles will coexist with a bulk gas 
phase only if their areas Ai = Af , with 

They then contain Nf B soluble gas molecules with 

«f = £^f^V (S) 



kT \7T/ \ 0" 



3 Nondilute Foams 



Now we consider "nondilute foams" in which bubbles press on one another and are 
distorted into nonspherical shapes f25 ]. Following Princen |21^p3|, we consider 



the compression of a previously dilute foam under an osmotic pressure IT. 



In a nondilute foam a typical bubble's interface consists of gently curved faces 
which contact adjacent bubbles, and highly curved regions at the Plateau bor- 
ders. Although bubble faces press on one another with a disjoining pressure, this 
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typically results in a negligible direct contribution to the free energy ||26|| ; we may 
assume that the disjoining forces only indirectly affects a bubble's free energy, by 
distorting its shape and increasing surface area. Put differently, throughout the 
foam the surface tension a is taken constant, independent of the volume fraction 
of gas present. 

Following Equation |] (the ideal gas law) we obtain A/i; = /tf — n s b for nondilute 
bubbles as A/tj = kT In (P//P 6 S ), where fi s b is the chemical potential of a bulk gas 
of soluble molecules subject to pressure P b . Were such a bulk gas to arise by 
coarsening, P b s would balance both the atmospheric pressure P and the osmotic 
pressure II; defining P = P + II we require P b = P (see figure [L]) p7[ . So for 
nondilute foam bubbles we have 

AiH = kTki(5L) (9) 



P 

As we compress the system with a semipermeable membrane, the previously 
spherical bubbles will distort in shape and the continuous phase will flow 
so that the additional pressure is evenly distributed amongst bubbles. (This 
contrasts with granular materials for example [^,|3(|.) In the simplest scenario of 
2D, monodisperse foams, the bubbles are compressed into monodisperse hexagons 
with 'rounded' corners, and equal internal pressures. However, compression of a 
general polydisperse foam will result in bubbles deforming into various shapes; a 
given bubble will have a pressure which depends not only on its volume (as for 
spherical bubbles) , but on the arrangement and pressures of all of its neighbours. 



The above discussion may be clarified by noting that bubble interfaces which 
press on one another, do so with a radius of curvature at most of order 1/Ri 
(where Ri is the radius of a bubble with the same volume in an uncompressed 
state). Hence pressure differences between bubbles are of order a/Ri. So if n 
exceeds a/Ri then the increase in bubbles' pressures will (to within terms of order 
a/Ri), be homogeneous throughout the foam pll . So we define the increase in a 
bubble's pressure above that of bulk gas (P b = P = P + n), by 

Ff = Pi - p - n (io) 

and expect Pf to be of order a/Ri. Here Pf reduces to the Laplace pressure 
of the ith dilute, uncompressed (n = 0), spherical 3D bubble, for which the 
increase in pressure above a coexisting bulk gas is Pf = Pi — P = 2a / Ri. But 
in a compressed state Pf is no longer the Laplace pressure, for the latter is cr/Vj, 
with Ti the radius of curvature at a Plateau border of droplet i. 

Since we consider an ideal mixture of soluble and trapped gases, then Eq. [TO] 
requires P* = P + H + Pf — Pj , which substituting into Eq. | gives 

/ P g_ p t\ 
A/^ = fcTm 1+ ~ I (11) 
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Figure 1: The pressure of a bulk gas P§ is negligibly increased by interfacial 
tension, so P fe s will balance (and hence equal) P + EL If the bulk gas arises by 
coarsening, so that one bubble grows to macroscopic size, its own trapped species 
contribute negligibly to the pressure. 

So since a bulk gas formed by coarsening has A/Zj = 0, a bubble can coexist with 
a coarsened bulk gas phase when P^ = Pj '. 

We note P? as defined by Eq. |10| is generally not easy to calculate from geometric 
considerations. Nonetheless the origin of Pf for a spherical bubble is geometrical, 
and Pf of a nondilute foam bubble is determined by packing geometry, so we 
refer to Pf as a bubble's "geometric pressure". To confirm the reasonableness 
of our arguments for the magnitude of Pp, the following section considers a 2D 
model for which Pj, II, and Pf are exactly calculable. 

We first clarify what it means for a foam to be "incompressible" . If Pf <C P + II 
then variations in a bubble's pressure are negligible compared with its actual 
internal pressure. So a bubble's gas density is approximately unaffected by its 
geometric pressure, and hence we may treat such bubbles as effectively incom- 
pressible. We emphasise that bubbles with Pf <C P + II may only be treated as 
incompressible with respect to coarsening (which changes their geometric pres- 
sures); such systems are not incompressible under changes in P or II (which 
changes their Laplace pressures as well). 
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4 A Simple 2D Model 



We now study a mono disperse, incompressible 2-dimensional foam, at osmotic 
pressure II. (In Appendix [A| an alternative argument confirms the results for 
an equivalent, but compressible model.) For II 7^ monodisperse bubbles will 
form an hexagonal array, with bubbles distorted into approximately hexagonal 
shapes but with rounded corners (see figure |2|). For simplicity the film thickness 
between adjacent bubbles is taken to be negligible. We note that an equivalent 




r 



Figure 2: A 2D monodisperse foam. Since bubbles are equivalent they have the 
same pressures and hence radii of curvature at their Plateau borders. 

system of osmotically compressed, monodisperse cylindrical emulsion droplets 
was considered by Princen [^I|-|23[| . Princen's calculation for the osmotic pressure 
is used later. 

Since the bubbles are taken to be incompressible, the osmotic pressure II does 
work by removal of the continuous liquid phase from bubbles' Plateau borders. 
Writing the area of liquid associated with each bubble as An (with An given as 
the sum of one third of the volume of liquid at each of its Plateau borders), then 
since we consider a monodisperse system with bubbles of area A4, the osmotic 
pressure is given by ]TB|, 
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where Li is the interfacial length of bubble i. (We keep a separate label i for each 
bubble, although they are identical, for clarity later on.) 



In Eq. |10]we defined Pf = Pi — P — XI. Since the monodisperse system considered 
here has hexagonal symmetry, by analogy with the Laplace pressure of a circular 
bubble we propose that P t G may be calculated by considering the increase in a 
bubble's interfacial length with an isotropic expansion at fixed liquid area An. 
So for this system we postulate that 



with (dCi/dAi)j\, H calculated for an isotropic expansion, and with Pf the same 
for all bubbles (since the system is monodisperse). Pf will be different for differ- 
ent lattice arrangements, and also varies (along with bubble shape), with liquid 
content. Only in the absence of bubble-bubble contacts will Pf equal the Laplace 
pressure, with Pf = a/Ri = a/ri. 

Since Pf is calculated at a fixed volume of liquid per gas bubble, we have fixed 
radius of curvature at the Plateau borders. Hence Eq. [13] gives 

where I is the length of the flat bubble-bubble faces (all equal). 

In Appendix |B] we obtain the exact expressions for the interfacial length, and 
area of a nearly hexagonal bubble as 



d = 61 + 2nr (15) 
Ai = T\r 2 + 6/r + ^5/ 2 (16) 



After differentiation and some algebra, Eqs [14], [y| and [16] give 

„n 6cr 



/36r 2 + 6^04* - vrr 2 ) 
We note that P G may also be written as 



(17) 



V Ai + A u 

where An = (2\/3 — n) r 2 is the area of liquid in Plateau borders, per bubble. 
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Princen fl21|,P2]| calculated the osmotic pressure of mono disperse, cylindrical emul- 
sion droplets by equating the work done by the osmotic pressure II with the in- 
crease in interfacial energy as droplets distort (at fixed droplet volume). For this 
quasi-2D geometry he obtained 



n 



a 
R, 




1/2 



1/2 



1-0 



(19) 



with the volume fraction of emulsion droplets, 0o their volume fraction at first 



contact, and Ri the corresponding radius. Princen's 2"2| calculation applies 
equally to our incompressible, monodisperse 2D foam: the area A{ and volume 
Vi of monodisperse cylindrical drops of length A obey Ai = £jA, Vi = AiA. 
Therefore O = ( 7r -R« 2 /2v / 3-R l 2 ) and = Ai/(Ai + An). In the variables used in 
this paper we have 




1/2 



'2^3\ 



7T 



\ 1/2 / Al \ 1/2 

J \Ai + AiJ 



2v/3 



V A + Ai 



-Ri 



and 




' 2x/3 
Ai + Ai, 



where comparison of Eq. |2(] with Eq. 
Eqs. M and EE Eq. HI becomes 



n 



1/2 



R; 



'2V3- 



7T 



1/2 



2v/3 / 
llj shows that (0/0 o ) 1/2 



Ai, 



1/2 



P^Ri- So using 



(20) 



Ri 

r 

(21) 



(7 

r 



P 



G* 



(22) 



So since Pj = P + er/r, we have the exact result that Pj = P + II + Pp, which 
for the system studied confirms both our physical argument for the pressure in 
an osmotically compressed bubble, and our hypothesis for this 2D monodisperse 
system that the geometric pressure should be calculated for an isotropic expansion 
at fixed Plateau border radius r (Eq. pf). 

From Eq. |17], in the dilute limit of circular bubbles (r 2 = Ai/rr) Pf — > oypHj yCAi, 



and in the dry limit of hexagonal bubbles P^ — ► cry 2y/3/ \fA~i. So since y 2\/3/ ^fix ■ 
1.05, we find the surprising result that Pf (and hence Pi), is only weakly affected 



by IT. We return to this in Section 



The area Af at which bubbles in a hexagonal packing may coexist with a bulk 
gas (so that P? = P + IT), obeys Pf(Af) = Pf(Af). So using Eq. [TJ| and 
Pj = N T kT/Ai we obtain 



A 



^Hi 



1 + 



1 + 4An/A% 



(23) 
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where A% = (1/ 'A^>){N T kT '/ 'a) 2 , is the with which an entirely dry foam 
with hexagonal bubbles (An = 0) may coexist with a bulk gas. We then obtain 
the number of soluble species at coexistence from Nf B = (P/kT)A B . 



5 A Condition for Stability 

Having established a general definition of geometric pressure (Eq. |T0|) , and shown 
that it corresponds (in at least one special case) to an isotropic expansion at fixed 
liquid content, we now use conservation of the total number of gas molecules and 
total number of bubbles to derive a criterion to ensure the formation of a stable 
distribution of foam bubbles. We take the bubble size distribution to be composed 
of two parts, a 'coarsening' part of the distribution, and a 'stable' part consisting 
of bubbles which have shrunk to a stable size at which they coexist with the 
coarsening bubbles. For a sufficient quantity of trapped species the assumption 
of the bubble distribution having a coarsening part is found to be inconsistent, 
enabling us to derive a stability criterion for foams. These arguments have strong 
similarities to those in ||, but are considerably generalised. 

We take raf , and as the number densities of bubbles overall, in the stable 
part of the distribution, and in the coarsening part of the distribution respectively. 
Conservation of bubble number (no coalescence) gives 

n o = n s + n c (24) 

where we note that nf and n% may be time-dependent. Nf B , the number of 
soluble species present in bubble i when coexisting with a bulk gas phase, is 
determined by P^ = Pj . However since relations between P^ and Vi may vary 
during coarsening (due to changes in the bubble's environment), we take Nf B as 
time-dependent. We define the following 

• N s0 = (Nf(t))i : The average number of soluble molecules per bubble. This 
is time independent. 

• N sB (t) : The average number of soluble molecules per bubble in the stable 
distribution. As coarsening proceeds, n^/n® — > and the bubbles in the 
stable distribution have iV/ — > N^ B . Then N sB (t) tends toward the average 
number of soluble molecules per bubble at coexistence with a bulk gas 
phase. 

• N sC (t) : The average number of soluble molecules per bubble, within the 
coarsening distribution, at time t. 
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As a coarsening foam tends towards an equilibrium state, conservation of the 
number of soluble molecules requires 



n° b N s0 = n s b N sB (t) + n^N sC {t) (25) 
Combining Eqs. [23] and gives 

n° b [N s0 - N sB (t)) = n% (N sC (t) - N sB (t)) (26) 



So if 

N sB {t) > N s0 (27) 

then Eq. |26| requires > n b (N sC (t) — N sB (t)). But since the larger, coarsening 
bubbles have N sC (t) > N sB (t), then n b = 0, and the foam must be stable against 
coarsening. Since N sB (t) depends on P,p(t), whose relation to V$ is not known, 
the derivation of a stability condition in closed form is not always possible, unlike 
the dilute case M. However, the above condition enables us to investigate the 



requirements for stability. This is done in Sections 5.1 and |6.1| below. 



5.1 Exact Stability Conditions 

For both dilute foams and the model of a monodisperse 2D foam, the geometric 
pressures are calculable exactly, giving exact expressions for Nf B . These are: 
Eq. §; Eq. @; and N° B = (P/kT)A B with Af given by Eq. |2S|; for dilute 3D 
foams, dilute 2D foams, and the monodisperse 2D model respectively. Stability 
conditions are calculated by averaging over the relevant equation for Nf B (where 
appropriate), and ensuring that Eq. |27| is satisfied. 

For dilute 3D foams the requirement is N s0 < N sB , with Eq. |^ giving 

which resembles, but generalises a result in . Similarly for dilute 2D foams we 
require iV s0 < N sB , but with Eq. [| giving 

Finally, since the model of 2D foams considers monodisperse bubbles we merely 
require that A^ s0 < Nf B , with Eq. p3| giving 



^ = | A s = |^ 1+ VL±i«| (30) 
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Recall that = (l/4 v / 3)(A^ T fcT/a) 2 so that N- B depends on a as expected. 
The area of liquid associated with each bubble becomes negligible compared to 
Atfi as the foam becomes increasingly dry. Note also that the dry monodisperse 
hexagonal foam is, even without trapped species, already dynamically stable with 
respect to infinitesimal volume changes of a single cell but not with respect to geo- 
metric reorganisation so as to create five-sided and seven-sided cells. Such a foam 
is, however, unstable with respect to homogenous cell shrinkage throughout the 
system (with the excess gas forming a bulk coexisiting phase), whereas sufficient 
trapped species, as calculated above, will restore full thermodynamic stability in 
this sense. We accept that the distinction between stability and metastability 
becomes blurred, especially as the geometry of the foam becomes more complex. 



6 Polydisperse Foams 

The edges of bubbles in a reasonably dry 2D foam (II ^> a/Ri), meet with an 
angle approximately equal to 2n/3. If the edges met with an angle of exactly 
27r/3, then bubbles would need to have exactly equal radii of curvature at their 
Plateau borders, and hence equal bubble pressures Pi = P + a/r. Since the 
curvatures between adjacent bubbles remain of order I/R4 <C l/l, then as a 
foam becomes increasingly dry, bubbles' pressure differences become increasingly 
small compared with the mean bubble pressure (Pj 3> Pf ~ a/Ri). So provided 
a polydisperse foam is sufficiently dry, then Pf ~ a/Ri and P® <C U. 

In a sufficiently wet and polydisperse foam, very small bubbles might reside 
within the Plateau borders of larger bubbles without mechanically experienc- 
ing an osmotic pressure II. Such a bubble in a 3D foam will be spherical with 
Pi = P + 1a/Ri and hence Pf = 2a /R { - U. For such bubbles Pf need not 
necessarily be of order a/Ri. In fact even for the same bubble size distribution 
and the same osmotic pressure II, a foam's volume can depend on the arrange- 
ment of bubbles it contains (see figure However, given a reasonable osmotic 
compression II, then a reasonably narrow distribution of Nf will be sufficient 
to ensure that such bubbles occupy a negligible volume fraction, and hence may 
be neglected when calculating the stability condition Eq. |S7|. This is shown in 
Appendix |C[ 



6.1 Stability Requirements 

To address the stability condition Eq. |27|, we need to know the size at which 
shrunken bubbles in the stable distribution will coexist with a bulk phase of 
soluble gas. We will assume that for a given II the distribution of Nf is sufficiently 
narrow that we may neglect any tiny bubbles residing wholly within Plateau 
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Figure 3: At a given osmotic pressure IT, the volume occupied by a weakly com- 
pressed polydisperse foam can depend on the way in which bubbles are arranged. 
The two foams shown both contain the same bubbles but the different bubble 
arrangements result in their occupying different volumes. 

borders of larger ones (see above, and Appendix 0). Then for a shrunken bubble 
in T> dimensions of volume Vi, we define a quantity Tj such that 

p ? = r*-^ (3i) 

and expect Tj ~ 1. (Any coarsening or rearrangements will make Tj time depen- 
dent.) For example, our model 2D foam has (from Eq. 



At 



2 ^ / 3\/ 1 " a a ( 32 ) 
which is graphed in figure [§. In this example Tj varies monotonically between 



\pii for an entirely dilute foam (circular bubbles), and y2v / 3 for an entirely dry 
foam (hexagonal bubbles), but always remains of order 1. 

We may calculate an approximate value for the volumes of shrunken bubbles by 
taking Tj as independent of V{ (although changes in a bubble's environment mean 
that Tj may change with time). Using N^ B = (P/kT)V B , and obtaining by 
equating PP = T i aVi /V with = NfkT/Vi, gives 

P ( N T kT\ V/(p ~ 1] 
N! B = ^( ^) (33) 



kT \ aTj 
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Figure 4: Variation of Tj (vertical axis), from y/n to y2\/3, with X = I jr. 
Increasing X corresponds to the bubble becoming increasingly dry, and X — ► 
corresponds to the bubble becoming circular in shape. 

In Appendix [D] we use Eq. ^ to investigate the stability requirement Eq. |27| 
in more detail. By assuming that and Nf are uncorrelated, it is shown that 
unless the average value T may increase without bound, then a stability threshold 
does exist. (That is, there will always be some number of trapped molecules per 
bubble beyond which coarsening will be prevented.) Equation |32| shows that, for 
a monodisperse 2D foam, T is bounded above even in the dry limit (in contrast 
to the Laplace pressure). Also since (P^) ~ a/R (see appendix |A|), then for P 
to increase without bound, Tj and Ri will need to be correlated in a very specific 
way. On balance, all these arguments suggests that f remains bounded as a 
polydisperse foam evolves. Hence our results suggest it should always be possible 
to osmotically stabilise a polydisperse foam by adding enough trapped species. 
As a rough estimate, stability requires a sufficient pressure of trapped gas that 
P T > cr/R, which for a ~ 10 _1 Nm _1 , and R ~ 10~ 5 m requires P T ~ 10~ 4 Nm~ 2 , 
ie of order 0.1 of atmospheric pressure. 
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7 Coarsening of Foams: Qualitative Behaviour 



Now we consider the coarsening of incompletely stabilised, non-dilute foams. Pre- 
vious work on coarsening of dilute emulsions in |J also applies to dilute foams 
mutatis mutandis. Here we concentrate on non-dilute foams in which bubbles 
impinge on one another and are distorted from their otherwise spherical shape. 

As with dilute emulsions @, the trapped molecules prevent bubbles from entirely 
disappearing. The resulting foam morphology and coarsening kinetics will be 
determined by two main factors. The first is the 'excess volume fraction' of 
dispersed phase, defined as the total volume fraction of gas which ultimately will 
coexist with a stable ensemble of shrunken bubbles, in osmotic and mechanical 
equilibrium with it: this is the amount of gas actually available for coarsening. 
The second factor is the ease with which shrunken bubbles may rearrange to allow 
larger bubbles to coarsen. 

The excess volume fraction of disperse phase, determines a foam's expected late- 
stage morphology (see figure [5]): 

1. No excess volume fraction: The foam is stable. 

2. Very low excess volume fractions: Larger bubbles are surrounded by a 'sea' 
of shrunken bubbles. Competitive coarsening between larger bubbles re- 
quires a gas flux through the sea of smaller bubbles. 

3. Very high excess volume fraction: Larger bubbles are decorated by collec- 
tions of smaller bubbles at their vertices, with their faces impinging on other 
large bubbles. 

4. Intermediate volume fractions: Structures between the previous two ex- 
tremes. 

In what follows, we focus on the case of low but nonzero excess volume fraction. 
Here coarsening of larger bubbles will require rearrangements of the shrunken 
bubbles to prevent the build up of excess elastic strains (which will otherwise 
halt coarsening, as shown below). We envisage four scenarios: 

1. Inviscid rearrangements: Bubble rearrangements occur easily, and with neg- 
ligible dissipation of energy. 

2. Viscous rearrangements: Bubbles can rearrange, but resist doing so and 
hence slow the rate of coarsening. 

3. Negligible rearrangements (elastic medium): Bubbles grow within an effec- 
tively elastic medium, which may eventually arrest coarsening. 
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Figure 5: Late stage foam morphologies. For low excess volume fractions (left) 
grown bubbles are well separated by large numbers of shrunken bubbles, whereas 
at high excess volume fractions (right) the grown bubbles impinge on one another 
with shrunken bubbles decorating the grown bubbles' vertices. 

4. Elasto-plastic rearrangements: There is a maximum yield strain beyond 
which rearrangements allow flow of the shrunken-bubble sea, but below 
which it behaves as an elastic medium. 



We expect rearrangements to occur easily in a sufficiently wet foam, but rear- 
rangements in a very dry foam to occur rarely or not at all. So we expect the 
scenarios from 1 — ► 3, to become more applicable as foams become increasingly 
dry. Since both rearrangements and diffusion of disperse phase are required for 
coarsening to occur, coarsening will proceed with a rate determined by the slow- 
est process. This contrasts, for example, with phase separation in a binary fluid, 
where the coarsening rate is governed by the fastest process (with diffusive coars- 
ening at early times, viscous hydrodynamic coarsening at intermediate times, and 
inertial hydrodynamic coarsening at late times ||34j| ). In the present case rapid 
coarsening may initially be limited by viscous forces, then later as coarsening 
slows, viscous forces will become negligible and coarsening diffusion-limited. 

Note that a similar classification of kinetic regimes may in part be applicable 
to the coarsening dynamics of foams containing no trapped species. It would 
be appealing, perhaps, to study this case in detail first, before addressing the 
situation where trapped species are present. However, the latter case is actually 
a lot simpler, at least in the case of low excess volume fractions (figure [5], left) 
considered here. This is because the actively coarsening bubbles are effectively 
decoupled by a sea of passive shrunken bubbles; this simplification is absent at 
excess volume fractions approaching unity (which recovers the unstabilised case). 
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8 Mean-Field Model: Inviscid Rearrangement 



We consider a mean-field model for the coarsening of a small excess volume frac- 
tion, in which grown bubbles are already sufficiently large that they contain a neg- 
ligible quantity of trapped molecules, and have an approximately spherical shape 
(figure 0). Hence large bubbles have Pf ~ 2a /R4 and Pf ~ P + 2a /R4. We re- 
strict ourselves to incompressible foams (in the sense described in Section |3|), and 
take shrunken bubbles to have an approximately constant size Vb = 47ri?B 3 /3, 
with Rb <C Ri- An average pressure P s {r,t) of soluble gas in shrunken bubbles 
at distance r from a grown bubble's centre is obtained by coarse-graining over 
bubbles at r, and we assume P s (oo,t) is the same for all grown bubbles. We 
at first take rearrangements to be inviscid, so that bubble growth is determined 
by the rate at which soluble gas diffuses through bubble-bubble interfaces in the 
shrunken bubble 'sea'; this assumption is relaxed in Section P. 
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Figure 6: We consider a grown bubble with R 3> R , and shrunken bubbles at a 
distance r from its centre to have an average pressure of soluble gas P s (r,t). 

Consider the flux of gas from shrunken bubbles at radius r to adjacent bubbles at 
radius (r + R B ), see figure Following the approach of Von Neumann (see p5| ). 
we take the flux of gas between two bubbles as proportional to both the pressure 
difference of their soluble gas, and the surface area through which the gas may 
pass. Defining a flux velocity per unit pressure K, we obtain an average volume 
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flux of gas from bubbles at r to bubbles at (r + R B ) of 

Jy(r, t) = KA-nr 2 (P s (r, t) - P s (r + R B , t)) r 



(34) 



Since R^> R B then (P s (r,t) - P s (r + R B ,t))/R B ^ dP s (r,t)/dr, so solving for 
a steady state gas flux with V.J = 0, we obtain a radial flux 



M*.t) = K*, r #R* p '< s >-r l °° ,i > 



and a droplet growth rate 



dR 
~dl 



KR 



B 



R 



AP s (oo,t)-^j 



(35) 



(36) 



where we used P S (R) = P + 2a/R, and wrote AP s (oo, t) = P s (oo, t) — P. For an 
incompressible system all bubbles have a volume per gas molecule equal to that 
of a bulk gas, v g = kT/P. So with a little algebra we get 



dR 
~dt 



KPR 



B 



R 



2avg\ 
kTR J 



(37) 



where e = (P s (oo, t) - P) /P. 



Previously the thickness of liquid films between bubble faces was taken as zero. 
Now we take such films to have a small but finite thickness d. Then assuming 
the rate of flux through liquid films to be diffusion limited, we may calculate K 
in terms of the diffusion constant for dissolved gas molecules D, the volume per 



gas molecule v g , P = P + II, and d. This gives |T3 



K 



Dv g C(oo) 



Pd 



So the growth rate may then be written as 



dR 
It 



Rb 
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R 
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(38) 



(39) 



Note that for a dilute foam [13 the factor R B /d is absent from Eq. |39|. The 
increase in growth rate is due to the reduced volume fraction of liquid through 
which gas molecules actually need diffuse, a reduction of order d/R B . 



9 Rate-Limiting Mechanisms 

We continue to study the mean-field model of Section ||, but no longer require 
inviscid rearrangements. When bubble rearrangements are inviscid there is dissi- 
pation from the diffusion resistance of dissolved-molecules diffusing through liquid 
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films between bubbles, and if bubble rearrangements are viscous there is also dis- 
sipation as bubbles rearrange. By equating the rate of dissipation with the rate 
of decrease in free energy, in Sections |9T2| and |T3] we obtain the order of 
magnitude for a droplet's growth rate. The greatest source of dissipation limits 
the coarsening rate, so by comparing the dissipation rates we can estimate (Sec- 
tion p.3|) when each type of coarsening will occur. A similar approach applied to 
emulsion rheology (not coarsening), is found in 



As an example, we firstly consider coarsening in the traditional LSW fTT| , |P2f 
scenario of a vanishingly small volume fraction of bubbles in a liquid (without 
trapped molecules). Here dissipation is from the diffusion resistance of diffusing, 
dissolved gas molecules. 



9.1 LSW Coarsening of Bubbles 

In a steady state the radial flux of dissolved gas molecules, at a distance r from 
the centre of a bubble of radius R is 

J{r) = ~ 40 

v g r 1 

We write J(r) = c(r)u(r), with c(r) and u(r), the concentration and average 
radial velocity of dissolved gas molecules at r respectively. Since c(r) ~ C(oo), 
where C(oo) is the concentration of dissolved gas molecules adjacent to a bulk 
gas phase, then 

R 2 R 1 

U r 777 — \~ 41 

v g C(oo) r z 

The number of dissolved gas molecules within spheres of radii r and r + dr is 
of order r 2 drc(r) ~ r 2 drC (oo). The average dissipation rate per molecule at r 
is (u(r) 2 , where ( is the viscous drag coefficient on a molecule of disperse phase 
moving through the liquid. Hence the total dissipation arising from diffusion 
to a drop is of order r 2 drC (oc) (u(r) 2 . Using the Einstein relation |36j ( = 



kT/D, Eq. O, and integrating, we obtain the rate of dissipation due to dissolved 



molecules diffusing between bubbles TSr>, as 

kT 

TS D ~ , 2 R 2 R 3 (42) 

Coarsening occurs to reduce interfacial energy, with a rate of reduction in energy 
of order aRR. So equating aRR with Eq. [|2| and rearranging, we get 

d _R D V 2 C(oo)a 

dt kTR 2 1 ' 
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in agreement with the traditional analysis of LSW [ 11| , [12| for coarsening of dilute 
emulsion droplets (as opposed to foam bubbles as studied here). 



9.2 Inviscid Rearrangements Revisited 

We now apply the method to the mean-field model of Section || where inviscid 



rearrangements were assumed. The argument in Section |9.1| gives the average 
velocity with which dissolved gas molecules will diffuse through liquid films, in 
Eq. ^1], and the dissipation per molecule in liquid films remains (u(r) 2 . But now 
dissipation only occurs in the liquid films between bubble faces. To integrate only 
over the volume of such films, the spatial volume element r 2 dr is reduced by a 
factor of d/FLs- So the dissipation due to diffusion resistance TSd is 

R B DC(oo)vj 



and we obtain 



Rb\ f Dv g C(oo) \ / ov g \ 



dR 

~dt r \ d J \ R I \kTRj 



(45) 



in agreement with the mean- field calculation (Eq. |39]). This gives a growth law 
for the mean droplet size, R ~ t 1 / 3 . 



9.3 Viscous Rearrangements 

We now consider the rate of dissipation due to viscous stresses in the thin liquid 
films, as shrunken bubbles rearrange. Adjacent bubble faces are again taken to be 
separated by a distance d, determined by the disjoining pressure between bubble 
membranes. Since the excess volume fraction is small, we assume that bubble 
growth results in fluid flow that is approximately radial. Taking v(r) as the 
velocity of the shrunken-bubble fluid, then incompressibility of that fluid requires 
V.v(r) = 0, so that in the spherically symmetric case there is a radial velocity 

v(r) = (46) 

at distance r from the centre of a bubble of radius R, growing with velocity R. 

Consider now two adjacent shrunken bubbles at distances r and r + Rb respec- 
tively from the centre of the growing bubble. The differing velocities at r and 
r + Rb will mean that bubbles must rearrange, and slide past one another. For 
r 3> R B the relative velocity of the bubbles is of order (dv{r) / 8t)Rb- The shear 
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rate of the liquid between bubbles is of the order of this relative velocity divided 
by the film thickness d. Within the film the viscous stress is therefore 

R B dv(r) 

Vl{r) ~ V— — — (47) 

where r\ is the viscosity of the continuous liquid phase, and j(r) is the shear rate 
in a liquid film at r. 

The volume-averaged viscous dissipation is dominated by the contribution from 
within the films [^(J and, per unit volume, is of order d/ Rb times the dissipation 
rate (of order r/7(r) 2 ) within each film. So we obtain the total dissipation rate 
due to viscous bubble rearrangements, TSy, as 

Using Eq. ^6] and integrating, this gives 

TS V ~ ^vRR 2 (49) 
d 

As in the previous calculations (Sections |9J] and |9^), the rate of decrease in the 
surface free energy of the growing bubble is of order aRR. Equating aRR with 
Eq. ^9|, and rearranging gives 

ad , . 

R ~ -— (50) 

and hence a linear growth law, R ~ t 1 . 

It is interesting to ask why the above argument and the resulting linear growth law 
does not apply to the conventional LSW coarsening of dilute emulsion droplets 
in a structureless fluid continuum. In emulsions the volume of a molecule in a 
droplet is similar to that in the continuous phase, so when a dissolved molecule is 
incorporated into a drop, the increase in the droplet's volume equals the volume 
of liquid previously displaced by the molecule. Hence the only displacement of 
liquid is that already accounted for by the Stokes-Einstein drag on the molecule 
as it diffuses through the continuous liquid phase. On the other hand, since 
the volume per gas molecule in a bubble is much larger than its volume when 
dissolved in solution, the viscous dissipation may be relevant for coarsening of 
dilute foam bubbles in a structureless fluid. (In this case, arrival of gas molecules 
at the surface of a growing bubble causes a net fluid flow radially outward.) A 
simple order of magnitude estimate gives 

TS V ~ 7]RR 2 (51) 

and R ~ cr/77. This differs from the case where the surrounding medium is a 
'shrunken bubble' fluid by the absence of the enhancement factor Rs/d. 
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9.4 Viscous or Diffusion Limited Growth? 



In a coarsening dense foam, both rearrangements and diffusion are necessary 
for bubble growth, so the growth rate will be limited by the greatest source of 
dissipation. So when TSd/TSv 3> 1 coarsening is diffusion limited, and when 
TSd/TSv <C 1 coarsening is limited by viscous dissipation. Here we estimate the 
ratio for plausible foam parameters. 



Comparing Eqs. HI and HU1 we see that 



TS D ( d\ 2 kT , 



TS V \RbJ DC{oo)v 2 Q r] 



R 2 (52) 



In terms of a molar concentration Cm, the molar volume vug of gas in bubbles 
at pressure P, and the gas constant Rq, we have 



TS D ( dV R G T 



-R 2 (53) 



TS V \RbJ DC M {oo)v Mg r] 

At room temperature and atmospheric pressure we take typical values of RqT ~ 
10 3 J, D ~ lfr 9 m 2 s-\ C M (oo) ~ 10 2 m" 3 , v Mg ~ l(T 2 m 3 , and r) ~ 10- 3 N m - 2 s, 
(eg. for C0 2 gas bubbles in water f3"8|). Taking R ~ 10~ 6 m, gives 



^~10»fAV (54 ) 

TSy \R B J ^ ' 

So that for micron-sized foam bubbles, viscous dissipation will be observed for 
djRs smaller than 10~ 2,5 . Thus for sufficiently small bubbles, thin liquid films, 
and high liquid viscosity, viscous limited growth may be observed (giving R ~ t). 
However at room temperature and pressure we expect diffusion limited coarsening 
to be more common, and moreover this will always dominate once the radius 
R of the coarsening droplets becomes sufficiently large. Note that for foams 
VMg = RgT/P, so TSd/TSv is proportional to P 2 /kT: at low pressures and high 
temperatures the prospect of observing viscous limited coarsening is increased. 
For example if P ~ 10 3 Nm~ 2 (1(T 2 atmosphere) then TS D /TS V ~ 10(d/R B ) 2 , 
which for d/Rs ~ 10 _1 predicts viscous limited growth. 

For dilute foams TSd and TSy are given by Eqs. |2] and ^T] respectively, so that 

TSd RqT 



TS V DC M (oo)v 2 Mg r) 



R 2 (55) 



Hence at room temperature and pressure, taking RqT, D, Cm(oo), r], VMg, and R 
as above, TSd /TSy ~ 10 5 and coarsening will be diffusion limited J39|. Writing 
VMg = RgT/P we find that viscous limited growth requires P 2 /RgT < 10 2 Nm 
which at room temperature requires P < 10~ 3 atmosphere. 
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10 Negligible Rearrangements: Elastic Medium 



Suppose we no longer allow bubbles to rearrange, so that the sea of shrunken 
bubbles acts as an elastic medium. By estimating the increase in elastic energy 
as a growing bubble changes volume, we show that coarsening will halt, with the 
foam now 'elastically' (as opposed to 'osmotically') stabilised. 

The assumption of negligible rearrangement requires detailed explanation, since 
in any foam (containing trapped species or otherwise) local rearrangement follows 
inevitable when certain conditions are reached at the junctions between adjacent 
cells. However, it is well known that, in the absence of coarsening, a foam can 
exhibit a finite macroscopic yield strain, below which there may be occasional 
local rearrangements but these are insufficient to allow plastic deformation of the 
medium whole fll 



Our physical picture is that the shrunken-bubble sea, which is fully stabilised 
against coarsening by virtue of the trapped species that it contains, can then 
offer elastic resistance to the growth of an isolated large bubble, even though the 
latter is at lower chemical potential of gas than other, less large bubbles that have 
coarsened elsewhere in the system. 

The elastic coarsening inhibition mechanism that we advance is thus only pos- 
sible because of the stability against coarsening of the shrunken-bubble sea. If 
the bubbles in this sea were themselves members of the coarsening population, 
then the rearrangement conditions at the junctions between them would be con- 
stantly met within it, and there would be an incessant rearrangement of bubbles 
caused by their local volume changes. Such a medium could not exhibit a yield 
stress in any conventional sense, and would be unable to elastically stabilise the 
growth of isolated large bubbles. This remark is consistent with von Neumann's 
theorem, which (in two dimensions) establishes that in the absence of trapped 
species, coarsening can never cease for any structure other than a perfect hexag- 
onal lattice. The theorem does not apply in the presence of trapped species 



and so the elastic arrest mechanism considered below leads to no contradiction 
with it. 

For definiteness we consider an initial state comprised of spherical ready-grown 
bubbles in an elastically unstrained sea of shrunken bubbles (see figure Such a 
system could be formed by osmotically compressing a previously dilute, partially 
coarsened foam. 

Consider a layer of shrunken bubbles initially at distance r from the centre of a 
grown bubble of initial radius Rq (see figure |7|). Then growth of the larger bubble 
from radius Rq to radius R will require the given layer to move so as to leave a 
volume equal to the change in the larger bubble's volume, which is 4^(i? 3 — Rq). 
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Figure 7: In the absence of rearrangements, bubble growth results in an increase 
in the shrunken bubbles' interfacial area. 



So a layer with initial inner radius tq must move so that ro is increased to r, with 

r 3 - r 3 = R 3 - R 3 (56) 
Hence we obtain a bubble layer's new position r(ro) as 

r(r ) = (r 3 + R 3 -Rl) 1/3 (57) 

In the absence of rearrangements, the linear extension H of bubbles in a layer 
with initial inner radius r is 

H = Au(r + R B ) - Au(r ) (58) 

where Au(r) is the distance by which shrunken bubbles at r will move (see figure 
0), due to the larger bubbles growth. So using Au(r ) = r(r ) — r , Au(r + R B ) = 
r(r + R B ) — (r + R B ), and Eq. |57|, we expand in terms of R B /r , which to lowest 
order gives 



r, 



H ~ R B -£ - 1 (59) 



r 



If we again use Eq. [57], then for small strains or large distances r from a bubble 
we may expand in terms of (i? 3 — Rq)/^, to obtain 



2„ fR 3 -Rl\ ^(R 3 -R 3 " 2 



H ~ --i2 fl — 5-31 +0 -^-31 (60) 
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Now consider the energy of extending or contracting along one axis of a single 
shrunken bubble, with its length changed from Rb to Rb + H . Restricting our- 
selves to small strains (or equivalently r ^> R), we may expand the associated 
increase in energy AEi(H) in terms of a power series in H, 

AE 1 = a ( ai H + a 2 H 2 + 0(H 3 )) (61) 

Since we consider an initial state with unstrained bubbles, equilibrium requires 
that both positive and negative values of H will give positive AE%, requiring 
ol\ = 0. So to lowest order AEi = a (a^-ff 2 ) ~ aH 2 . 

Initially we have of order rQdr /R B shrunken bubbles between spheres at r and 
r + dr . So the change in elastic energy due to a large bubble's growth or 
shrinking is 

< 62 > 

which upon substitution of Eq. |K| AE\ = aH 2 , and integrating gives 



DC 



AE ~ a k^- H (63) 

The asymptote of AE ~ R 6 for R ^> R means that in the absence of bubble 
rearrangements the elastic energy in the surrounding foam will ultimately prevent 
coarsening, and elastically stabilise the system. 



11 Finite Yield Strain 



The approximation of AEi ~ aH 2 is strictly only valid for large Tq or sufficiently 
small strains. For large strains we might expect plastic rearrangements to occur. 

Consider the following simple model: the shrunken bubbles may elastically sup- 
port a maximum (yield) strain y*, beyond which macroscopic rearrangements will 
occur. Then 

y* ~ H*/R B (64) 

with H* the plastic threshold for bubble rearrangement. In the absence of rear- 
rangements Eq. ISO gives 

I/? 3 - R 3 \ 

H ~ R B l 3 01 (65) 
r o 

so that H* implicitly defines a radius Tq within which rearrangements occur, but 
beyond which the medium behaves elastically. So from Eq. ^5], we may define 
r* = (R B /H*) (\R 3 - R 3 \), which since H* ~ y*R B gives 

rf - (66) 
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For a plastic region to exist around a bubble, we require Tq > R; Eq. |6| then 
requires \R 3 — R^\ > y*R 3 . 

Whether sufficient growth is possible for this threshold to be reached depends 
on the radius R' at which bubble growth would halt, in a purely elastic system 
as considered in the previous section. If the inequality \R' 3 — R^\ > y*R' 3 is not 
satisfied (as occurs for large enough y*) one requires a large fluctuation, either 
thermally or in the initial condition, to establish a plastic zone around a droplet, 
allowing it to coarsen further. In that case the foam is at least metastable and 
in practice this may be sufficient. 

To get a stability condition, we consider next the elastic energy stored around a 
droplet that has grown from size Rq to size R. Taking deformations within any 
plastic region to be of order H*, we may calculate the elastic energy as 

r ° r l dr o tt*2 . f°° r%dr 

~ +a /. ~w 

B Jr o U B 



ae -j:,^t h - 2 ^lw h{ ^ (67) 



giving 



a 



R 



B 



AE ~ y*^- ( R 3 - R 3 - y*R 3 ) + y*-- R 3 - R 3 (6 



a 



R 



B 



where, since y* < \R 3 — Rq \ / Rq, both terms on the right are positive. 

We now assume that a grown bubble adopts a state of mechanical equilibrium 
with the surrounding shrunken-bubble sea, on a timescale fast compared to any 
coarsening process. Then for a reversible fluctuation in bubble volume by dV, the 
work done by the bubble's gas equals the work done on the bubbles surrounding 
environment. Hence 



dV dV 



So differentiating Eq. with AE given by Eq. |53] for \R 3 — Rq\ < y*R 3 and by 
Eq. |68| for |^ 3 - _Rq| > y*R 3 , we get 



a 

p + n + — 

Ra 



R 3 -F& 



Rb ~ W) 



a \R 3 -R 3 \<y*R 3 

;; ^ R 3 - R o I C>3 _ D3| > *D3 



(70) 



This means that any bubble that has successfully coarsened to become a bulk gas 
phase (whose pressure clearly cannot be infinite as R —>■ oo) has overcome the 
plastic threshold and has a pressure obeying 

P b -P-U~y*^- (71) 
Kb 



Now we are able to estimate a stability threshold for a foam that combines a finite 
amount of trapped species with a finite yield strain (under the simplified initial 
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condition we chose above, that the initial state has a dilute array of partially 
grown bubbles in an unstrained, shrunken-bubble sea). In this case the lowest 
possible free energy obtainable by coarsening is that of a bulk gas at pressure P&. 
So unless the initial pressure Pj of grown bubbles exceeds Pb, then the foam must 
be stable . An initially unstrained bubble has 

P . - P - n ~ ^- (72) 

Kq 

so instability requires Pi > Pb, and hence 

n Kb 

So since Rb -C Rq, such a foam may only be unstable if y* <C 1. For example, 
considering foams formed by the osmotic compression of a dilute partially coars- 
ened foam, then if R is small enough that Rq < Rb/u* then coarsening may 
occur, but if it has already coarsened sufficiently that R > Rb/v*, then after 
osmotic compression it will be elastically stabilised. 



12 Conclusions 



We have considered dilute foams with spherical bubbles, and nondilute foams in 
which bubbles are compressed by an osmotic pressure II (which distorts bubbles 
and increases their gas pressure). Prior to osmotic compression, variations in 
bubble pressures are of order a/Ri, so bubbles may only support large osmotic 
pressures (II ^> a/Ri) by deforming and pressing on one another. Pressure 
differences between such osmotically compressed bubbles generally remain limited 
to a/Ri (with a possible exception in the case of extreme polydispersity in initial 
droplet size). If coarsening occurs, a coexisting bulk gas is created that has 
pressure P& = P + II, with P the ambient atmospheric pressure. We defined the 
excess in a bubble's pressure above that of such a bulk gas by Pf = Pi — P — II, 
and argued that Pf ~ a/Ri. This was confirmed quantitatively and explicitly for 
a monodisperse 2D model of foam in Section || P^ , which we call the geometric 
pressure, provides the driving force behind coarsening in the case of a compressed 
foam or emulsion; in general it is very different from the Laplace pressure but 
reduces to it in the dilute limit. For dilute 2D and 3D systems and also for 
a concentrated but monodisperse 2D foam, Pp can be found explicitly. This 
allows one to find rigorous conditions (Eqs. [2^, and |3D|) for their stability. 
More generally we have argued that for a reasonably large II and a reasonably 
narrow distribution of Nf, osmotic stabilisation should generally be possible, and 
typically requires Pj > a/Ri. 

The morphology of an insufficiently stabilised foam (containing inadequate levels 
of a trapped species), is determined by the 'excess volume fraction' of dispersed 
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phase; this will form a bulk gas phase if the system is allowed to reach full equi- 
librium (i.e., subject only to the constraint that the trapped species stay in their 
designated droplets and coalescence is absent, which with trapped species present 
will ensure conservation of droplet number). The kinetics of coarsening will be 
determined either by the rate of diffusion of gas between bubbles or the rate at 
which shrunken bubbles rearrange, with the slowest process determining the bub- 
ble growth rate. We studied coarsening of a small excess volume fraction, finding 
diffusion limited growth to result in R ~ t 1 / 3 (as opposed to R ~ t l l 2 as arises in 
foams without trapped molecules [j35]), and viscous limited growth to give R ~ t. 
For typical parameters we expect diffusion limited growth to be more common, 
with viscous limited growth more likely to be observed if the liquid viscosity is 



high, the temperature is high, and the pressure is low [fill . Under conditions of 
negligible bubble rearrangements, the buildup of elastic stresses among the sea 
of shrunken, stabilised bubbles will arrest coarsening of well-separated large bub- 
bles and 'elastically stabilise' the entire foam. This can sometimes occur when 
there is a finite yield strain; whether a foam will be elastically stabilised is then 
determined by the foam's initial state (and the yield strain y*). Thus osmotic 
compression of a partially coarsened dilute foam, may halt coarsening and elas- 
tically stabilise the foam; the further the foam has coarsened initially, the more 
likely this is to occur. 

Our work applies, substantially unchanged, to concentrated emulsions. As with 
nondilute foams the disjoining energy between repelling droplets in emulsions 
may usually be neglected PBJ, and the osmotic pressure of trapped molecules 



in droplets equals the partial pressure of trapped gas molecules in bubbles (so 
long as both are treated as ideal ||[TB|]). Hence if we replace the volume per 



gas molecule v g = P/kT with the molecular volume of disperse phase in droplets 
v b , and also replace Nf B = {kTjP)V B with Nf B = (l/v b )V B , then most of the 
above work applies equally well to concentrated emulsions. 



A Force Balance 



We adopt a similar approach to one previously used by Princen pl|,p2|, which 
will enable us to determine the order of magnitude of Pf when bubble sizes are 
polydisperse. Unlike the one in the main text, this approach is equally valid for 
compressible and incompressible bubbles. 

Consider a semipermeable membrane with a shape which follows the top surface 
of a line of hexagonal bubbles (see figure §), so that bubbles remain hexagonal in 
shape. Then since the membrane doesn't move, the total force on the membrane 
due to II + P must balance the total force on the membrane acting from below. 
Then we consider a typical bubble adjacent to the membrane (see figure §). Such 
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Figure 8: Imagine a flexible semipermeable membrane initially between bubble 
layers, which subsequently hardens and enables you to remove bubbles above the 
membrane, while maintaining a pressure II on the bubbles below. At equilibrium 
the forces above and below the membrane must balance. The same thought 
experiment may also be applied to polydisperse bubbles. 

a bubble will have a geometry as in figure || but with the bubble rotated by 
7i / 6 to leave two flat edges vertical. The force acting from below the membrane 
is Pi multiplied by the projection onto the horizontal axis of the flat bubble 
edges(p x 2(i\/3/2)), plus P multiplied by the projection of the Plateau borders 
onto the horizontal axis (P x 2r). The force acting from above the membrane 
is simply P + II multiplied by the bubble's width 2(/\/3/2) + 2r. Then for the 
forces above and below to balance we require [H2[ 



P,2^ + P2r = (II + P) + 2r) (74) 

Writing Pj = P + a/r and rearranging, we obtain 

n = - (75) 

which writing in terms of A and A\ and comparing with Eq. [18| gives II = a jr—Pf 
as before. Note that since the method may be applied to compressible systems, 
both Pj and Pf for the monodisperse 2D model will remain the same for both 
compressible and incompressible bubbles. 

We now generalise the argument to a polydisperse system. The projections of 
the bubble faces and plateau borders onto the horizontal axis are written as lp 
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requires 



and If respectively. Then considering a suitably shaped membrane, force balance 



{Pii-% + P(if)i = (n + p) + (ifh) (76) 

So writing Pj = P+a/ri and restricting ourselves to fairly dry foams with ~ rj, 
then (lf/ri}i ~ (lf)i/r where r = (rj)j, and we can obtain 

n^*-* <%- (77) 

Now we note that since if ~ r and if* ~ P, where P = (vC4/7r)j, then II ~ 
a/r - a/R so that if P? = P< - P - II then (Pf ) ~ a/P. 

In 3D the same argument but with projected (Af) ~Pxr, and (Af) ~ P 2 
leads to II ~ a/r — a/R which for Pf = Pi — P — II gives 

(If) ~ a/P (78) 

Finally we note that in higher dimensions the equivalent argument will continue 
to give (PP) ~ cr/P. 



B Geometry of Monodisperse, 2D Bubbles 



The area of a hexagon of side I is six times the area of the equilateral triangles of 
which it is composed. The area of each triangle is (l/2)h, with h = Zsin(7r/3) = 
ly/3/2. Hence we obtain the area of a hexagon as A = (3 v / 3/2)/ 2 . Its interfacial 
length is simply given by C = 61. 

We now describe a nearly hexagonal bubble in terms of the length of the flattened 
faces I, and the radius of curvature of the plateau border r (see figure ||). Then 
the area A of a nearly hexagonal bubble is given by the sum of six sections of a 
circle, each of angle tt/3 and radius r, plus six rectangles of length I and height 
r, plus the area of a hexagon of side length I. So 

A = Trr 2 + Qlr + ^-l 2 (79) 
and the surface length obeys 

C = 61 + 2vrr (80) 



C Bubbles in Plateau Borders 



What is required for states such as that in figure |3| to exist? Using Pf = P + 
2a I ' Ri — Pj T and the ideal gas law (for the trapped gas), the pressure of soluble 
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Figure 9: Nearly hexagonal foam bubble with flat edges of length 1, and radii of 
curvature r. The perpendicular distance between two flat edges is + 2r. 



gas in a spherical bubble residing in a Plateau border is P? = P + 2a /R4 — 
Nf kT/(4ir/3)Rf, which has a maximum value of P + 4(87r) 1 /V 3 / 2 /9(A;TiVf ) 1 / 2 . 
So in a system which coarsens with P s — > P+IT, the existence of spherical bubbles 
in Plateau borders requires 




lira 3 



>P + n (81) 



kTNf 

Hence such bubbles will need a sufficiently small quantity of trapped species that 

rp /8\ 2 2a 3 1 

N ' £ (5) ww (82) 

which writing IT in terms of the Laplace pressure of a bulk gas, II = 2a /r&, 
requires 

T /8\ 2 ar 2 err 2 



We may compare Nf with that for a typical stable (and compressed) bub- 
ble of radius R, containing N T trapped molecules. Then since P T (N T ,R) = 
N T kT/(47r/3)R 3 ~ a/R one requires 

RT ~JF (84) 

A spherical bubble in a Plateau border has 

NT < (85) 

Considering II for which r^/R ~ 10 _1 or smaller, then we require < 10 _2 A^ T . 
So for a reasonable II of order a"/r&, and a reasonably narrow distribution for 
A^, then there will be a negligible volume fraction of bubbles residing in Plateau 
borders. 
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D Stability Requirements 



We examine the form of N^ B (t), the number of soluble species per bubble below 
which a stable foam is ensured. Noting that N* B (t) = (N^ B {t))i, we write it out 
in full to obtain, 

P (kT\ V/{V ~ 1] //JV r \ B/M \ 
<^)>* = ^ - ( (¥\ ) (86) 



kT \ a J 

We note that V > 2 so that 



(0 



— > v/{v-i) (87) 



we assume that variations in P- G are independent of iV? . Then using Eqs. BB] and 



We assume that I\ and Nf are uncorrelated, that is for a given bubble volume 



37], we have 

6 , hT \V/{v-x) ( T v/(v-i), 

<™< s If (f ) 1 |^r < 88 > 

In Section [A] we found that (P G ) ~ &/R, which suggests T is of order a constant 
To Then since 

P fkT\ V/{V - l) (N TV/(p - 1) \- 

W B ^' W (t) i %^r i < 89 > 

a foam formed with an initial number of soluble species N s0 such that 

P fkT\ V/{V ~ 1] (NTW- 1 )) 

— ( — ) ^ h > iv s0 r go) 

kT { a ) fV/iv-i) ^ ^) 

will be stable. Similarly if T is bounded by T max , then a stability condition is 
obtained by replacing f c with f max in Eq. |90|. Alternately, if f is a decreasing 
function of time, then f (0) > T(t) and hence 1/f {t) v ^ v -^ > l/T(0) v / ( - v ~ 1 \ so 
provided 



P (kT 



v/(v-i) lATT v/{v-x\ 

> N s0 (91) 



kT \ a) T(0) v /( v -V 

then we could again guarantee a stable foam. In contrast to the above cases, 
if f increases without bound, then N sB (t) will decrease without bound, and a 
stability condition will never strictly exist. The only way a stability condition 
can fail to exist is if f can increase in time without bound. 
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E Table of Notation 



Ai, An, Af , area of two-dimensional bubble i; area of surrounding liquid 

(within hexagonal unit cell); area at which such a bubble (containing trapped 
species) coexists osmotically with a bulk gas phase; the same, but for a completely 
dry hexagonal foam. 

C(R): equilibrium concentration of dissolved gas adjacent to a bubble of radius 
R. 

c(r) : mean concentration of dissolved gas within the continuous phase at distance 
r from the centre of a bubble. 

D: diffusivity of dissolved gas within the continuous phase. 
V: dimensionality of space. 

d: thickness of films of continuous phase separating faces of adjacent bubbles. 

H, H*: local extensional or compressive strain of bubbles within a sea of shrunken 
bubbles; its maximum value before yielding occurs. 

Jy: volume flux of soluble gas around coarsening bubble. 

K: flux/pressure coefficient for gas permeation across a fluid film. 

Li. the perimeter length of bubble % in a two-dimensional foam. 

/: length of a flat interface between adjacent bubbles in a two-dimensional hexag- 
onal monodisperse foam. 

Nf, Nf, N* B : number of trapped molecules in ith bubble; number of soluble 
molecules within ith bubble; the number of soluble molecules within ith bubble 
when coexisting with a bulk gas phase. 

N s0 , N sB , N sC : average number of soluble gas molecules per bubble; average num- 
ber of soluble molecules per bubble within the stable (noncoarsening) population; 
average number of soluble molecules per bubble within the (unstable) coarsening 
population. 

n°,nf,n^: total number density of bubbles in the system; number density of 
shrunken bubbles; number density of coarsening bubbles. 

P: ambient pressure at which a foam is maintained. 

P: defined as P = P + IT. 
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P 6 S : partial pressure of soluble gas in a bulk gas phase created by a coarsening 
process, and equal to the total gas pressure in this phase. 

Pi, Pf, Pf: gas pressure within ith bubble; partial pressure of soluble species 
within it; partial pressure of trapped species within it. 

Pj G : geometric pressure of ith bubble (see text for definition of geometric pres- 
sure) . 

P s (r, t): mean partial pressure of soluble gas in shrunken bubbles at radius r from 
a growing bubble. 

Ri, R, R , R(t), R', Rb- radius of ith bubble; mean bubble radius; initial radius of 
a coarsening bubble; its radius at time t; radius at which its growth is arrested 
in absence of rearrangement; radius of a shrunken bubble at or near coexistence 
with a bulk gas phase. 

R G : the gas constant. 

r, tcTq: radial distance of a given layer of shrunken bubbles from the centre of 
a growing bubble; its initial value; the latter quantity for bubbles whose local 
strain coincides with the yield strain of the shrunken bubble sea. 

Ti. radius of curvature at Plateau border of ith bubble. 

Sd,Sv' entropy production due to dissipative and viscous fluxes during coarsen- 
ing. 

u(r),v(r): radial velocity of dissolved gas molecules at distance r from the centre 
of a coarsening bubble; radial velocity of shrunken-bubble fluid. 

VijVj 3 : volume of bubble i; its volume at osmotic coexistence with a bulk gas 
phase. 

v g ,VMg' molecular volume of disperse phase species; its molar volume. 
y*\ yield strain of the shrunken bubble sea. 
cti, a.2- numerical coefficients (dimensionless). 

AE, AEi(H): change in elastic energy of shrunken bubble sea on growth of a 
coarsening bubble; change in energy of a single shrunken bubble due to deforma- 
tion H. 

A/if chemical potential difference, /if — /if. 

e: dimensionless supersaturation (P s (oo,t) — P)/P. 
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rj: viscosity of continuous phase. 

0, 0o : volume fraction of bubbles; its value at first contact. 
Ti,T: for the ith bubble, the value of P^V^ V a\ its mean value. 
7: shear rate. 

UliHl'. chemical potential of soluble disperse phase in ith bubble; chemical po- 
tential of a bulk phase of the same species (at pressure P). 

II: osmotic pressure in a foam. 

a: surface tension of films in a foam. 
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